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FLIGHT INVESTIGATION AND ANALYSIS 
OF THE WING DEFORMATIONS ON A SWEPT-WING BOMBER 
DURING ROLLING MANEUVERS 
By Alton P. Mayo and John F. Ward 

SUMMARY 


The results of deflection measurements made at six spanwise stations 
on the wing of a swept-ving jet bomber (Boeing B-47A) during rolling 
maneuvers are presented in the form of coefficients expressing the deflec- 
tions due to the aileron loads, the sideslip loads, the wing-flapping 
inertia loads, and the rolling-velocity loads on the airplane. The pro- 
cedures used to obtain the coefficients are presented along with compari- 
sons of the experimental deflections and twists with those obtained from 
theoretical calcinations. Comparisons of measured and computed deflec- 
tions using available methods indicate good agreement. Detailed explana- 
tions of the least-squares procedures used in the analysis of the flight 
data and the methods used in the theoretical calculations are given in 
the appendixes. 


INTRODUCTION 


The design of a high-aspect-ratio swept-wing structure requires not 
only the application of advanced methods of stress analysis and the deter- 
mination of the wing stiffness, but also the determination of aerodynamic 
and inertia spanwise load distributions that would occur under specified 
flight conditions. In each field the computations are not only lengthy 
but are subject to a sufficient number of assumptions so that the end 
result obtained when the results of the separate fields are combined are 
sometimes in doubt. 

In order to secure an indication on the accuracy obtainable in each 
field, as well as of the accuracy of the end results, an extensive 
research program was undertaken in which a Boeing B-4-7A was used to obtain 
the experimental data. 

In one phase of the research program, deflections of various surfaces 
were measured for the purpose of comparing the measured deflections with 
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computed values. Initial results on this phase have been given in ref- 
erences 1 and 2 which give the deflection influence coefficients for the 
wing and the analysis of the wing deflections obtained in symmetrical 
flight. 

The present report is concerned mainly with the analysis of wing 
deflections obtained under unsymmetrical flight conditions. The unsym- 
metrical deflections were analyzed to obtain wing deflection coefficients 
expressing the deflections due to rolling velocity, aileron deflection, 
wing flapping, and sideslip. Comparisons are made between the experi- 
mentally determined deflections and those calculated theoretically. The 
procedures by which the flight data were analyzed for the various effects 
and the methods used in the theoretical calculations are given in the 
appendixes . 


SYMBOLS 


b' 


N 


Djji 


wing span minus fuselage width, in. (fig. 3) 

normal load factor at the tip of the rigid wing due to a 
unit pitching acceleration 

normal load factor at the tip of the rigid wing due to a 
unit rolling acceleration 

number of equations in least-squares solution 

airplane normal load factor at airplane center of gravity, 
positive when inertia loads are downward (n = 1 in level 
flight) 

normal load factor measured at the wing tip, positive when 
inertia loads are downward 


■“Tf 


P 

P 


normal load factor at wing tip due to wing flapping, positive 
when inertia loads axe downward 


airplane ro lling angular velocity, positive for right roll, 
radians/sec 

airplane rolling angular acceleration, positive for increasing 
positive rolling velocity, radians /sec^ 

airplane rolling angular acceleration due to loads on horizon- 
tal and vertical tail, positive for increasing positive 
rolling velocity, radians/sec^ 
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<1 

q 

x 

y 


Y' 


z" 

z 


Znj3 



airplane pitching angular velocity, positive for airplane 
nose pitching up, radians/sec 

airplane pitching angular acceleration, positive for increasing 
positive pitching velocity, radians/sec^ 

streanwise distance from intersection of front span center line 
and airplane center line, positive rearward, in. (fig. 3) 

lateral distance from airplane center line, positive left, in. 
(fig. 3) 

lateral distance from airplane center line less one -half the 
fuselage width, positive left, in. (fig. 3) 

vertical dimension, (measured from and perpendicular to the 
top of the fuselage) positive downward, in. 

total optigraph target deflection measured from wing-drooped 
position (ground zero), positive upward, in. 


target deflection due to wing inertia per unit airplane normal 
load factor, positive upward, in./n 


target deflection per degree of sideslip at unit airplane 

normal load factor at the center of gravity, positive upward, 
in. /deg 


target deflection due to wing flapping, positive upward. 



Zo 




z 


q 


target deflection due to the wing airloads when the summation 
of the aerodynamic loads on the airplane is zero, positive 
upward, in. 

target deflection per unit airplane rolling velocity, positive 
upward, in./radian/sec 

target deflection per unit airplane rolling acceleration, 
positive upward, in./radian/sec^ 

target deflection per unit airplane pitching velocity, posi- 
tive upward, in./radian/sec 

target deflection per unit airplane pitching acceleration, 
positive upward, in./radian/sec^ 
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z roll 


target deflection due to combined unsymmetrical loadings 
associated with the rolling maneuver, positive upward, in. 


Z ' roll 


3 AIL 


P 

6 AIL 

a 


a ^AXL 


target deflection due to the combined unsymmetrical loadings 
plus the effect of the zero-lift loads and wing droop, 
positive upward, in. 

target deflection per unit aileron deflection positive upward, 
in. /deg 

angle of sideslip, positive for right wing forward, deg 

left wing aileron deflection, positive downward, deg 

wing section streamwise angle of attack, positive leading edge 
up, radians 

wing streamwise angle of twist due to unit aileron deflection, 
radian/deg 

wing streamwise angle of twist due to unit wing flapping load 
factor, radian/nj 


a wing streamwise angle of twist due to unit sideslip angle, 

radian/np 

a, wing streamwise angle of twist due to unit rolling velocity, 

^ radian/radian/sec 

Matrix symbols: 


L J 

{ I 


[] 



row matrix 
column matrix 
square matrix 
rectangular matrix 
inverse matrix 
diagonal matrix 
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AIRPLANE AND TESTS 


5 


The airplane used in the test was a Boeing B-47A. (See figs. 1 
and 2.) The c hang es in the test airplane configuration from the standard 
airplane were the installation of (l) an airspeed measuring boom and 
fairing on the nose and ( 2 ) an external canopy, housing the deflection- 
recording instruments, mounted atop the fuselage approximately at the 
intersection of the airplane center line and the wing 38 -percent -chord 
line. 


The flight-test data used in this paper pertain to four aileron-roll 
maneuvers flown during the B-VfA flight research program conducted at the 
NACA High-Speed Flight Station at Edwards, Calif. 

The four aileron-roll maneuvers consisted of left and right rolls 
at Mach numbers of 0.60 and O. 7 I. The specific values of aircraft weight, 
center -of -gravity position, Mach number, altitude, and dynamic pressure 
axe included in table 1 for each roll maneuver. 


INSTRUMENTATION AND ACCURACY 


The instrumentation pertinent to the results presented in this paper 
consisted of a pitch turnmeter, roll tummeter, yaw tummeter, altimeter, 
airspeed recorder, an accelerometer at the center of gravity of the air- 
plane, an accelerometer at the left wing tip, an aileron control-position 
recorder, a sideslip angle recorder, and an optigraph system for recording 
the wing deflections. 

Corrections were made to the recorded center-of -gravity-accelerometer 
data for the small displacement from the airplane center of gravity. 
Corrections were made to the sideslip data for the effects of yawing veloc- 
ity and induced flow. The induced-flow effects were small, that is, on 
the order of 5 percent, and were based on estimates of the flow around the 
fuselage and boom. All instruments used were of the standard NACA photo- 
graphically recording type with the exception of the optigraph system which 
was designed especial l y for the complete test program on the airplane. 

The wing optigraph system consisted of eight target lamps on the left 
wing and optical recording instruments located atop the fuselage approxi- 
mately at the intersection of the 38 -per cent -chord line and the center line 
of the fuselage. (See fig. 3 .) To facilitate recording the deflections 
optically in the daylight, high-intensity infrared light sources were used 
in combination with infrared sensitive recording film. The optigraph sys- 
tem was calibrated through the use of a calibration stick, with 12 lamps 
on it at 6 -inch intervals, held vertically at each target station during 
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the calibration. All inflight measurements were made with reference to 
the wing -droop position with the airplane on the ground and with the 
outrigger gear clear. 

The instrument sensitivities, locations, and the estimated accura- 
cies of measurement are given in table 2 . 


METHOD AND RESULTS 


In a maneuver combining pitching, yawing, and rolling, the final 
load distribution may be considered as arising from the superposition 
of various types of load distributions. In the analysis of the flight 
data and in the theoretical calculations of this paper, the final load 
distribution is considered as having the eight components enumerated 
subsequently. The wing structural deflections are assumed to be line- 
arly related to these component load distributions . The following sym- 
metrical distributions are considered: 

(1) A zero lift aerodynamic loading at zero airplane load factor 
necessary to balance fuselage and tail loads. The shape of this distri- 
bution is affected by actual geometric wing twist as well as by twist due 
to the interference effects of the fuselage and nacelles. 

(2) An additional type of load distribution proportional to airplane 
load factor. This distribution has aerodynamic and inertia subcomponents. 

( 3 ) An airload distribution associated with pitching velocity which 
is a result of the wing -angle -of -at tack change caused by the velocity of 
the wing perpendicular to the airstream. 

(4) A loading associated with the angular acceleration in pitch - 
mainly an inertia loading. 

The following unsymmetrical distributions are considered: 

(5) A load distribution which is due to rolling velocity and which 
has aerodynamic and inertia subcomponents. The aerodynamic subcomponent 
is a result of the wing -angle -of -attack change caused by the velocity of 
the wing perpendicular to the airstream. The inertia subcomponent is a 
result of the airplane rolling acceleration associated with the rolling- 
velocity aerodynamic load. 

( 6 ) A load distribution which is due to aileron deflection and 
which has aerodynamic and inertia subcomponents . The inertia subcomponent 
is a result of the airplane rolling acceleration caused by the aileron 
airload. 
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(7) A load distribution which is due to sideslip and which also has 
aerodynamic and inertia components. The inertia component is a result of 
the airplane rolling acceleration caused by the sideslip airloads. 

(8) A load distribution which is due to wing flapping - mainly an 
inertia distribution; for this airplane the distribution has the frequency 
of the wing first unsyrametrical bending mode. 


Analysis of Flight Data 

In order to illustrate the type of flight <jLata used in the analysis, 
time histories of the deflection of target 9 and the associated airplane 
motions are shown in figures k, 5; 6> and 7 for the four roll maneuvers 
investigated. 

The procedure by which such measurements were reduced to coefficients 
expressing the deflection at a point due to the various types of loads is 
as follows: In the range where wing section lift and pitching moment due 

to aileron are linear, the wing deflections caused by the aileron 
aerodynamic -plus -inertia load component (item 6) may be expressed as 


Z = Zc 


3 AIL 


3 AIL 


( 1 ) 


Also, from consideration of the loads involved, the wing aerodynamic- 
plus-inertia deflections caused by the rolling velocity (item (5)) of the 
airplane may be expressed as 


Z = Zpp 


( 2 ) 


The deflection due to wing flapping (item (8)) is a result of an 
inertia effect and may be expressed as 




The load factor rv^ at the wing tip caused by wing flapping was obtained 

from the load factor measured at the wing tip by eliminating the effect of 
the airplane center-of -gravity normal load factor and the effect of. rolling 
and pitching acceleration. This is expressed as 
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= n-r - (n + K^p + K^q) 


00 


where the expression within parentheses represents the computed normal 
acceleration at the wing tip for the rigid airplane. Time histories of 
the load factor computed from equation (4) are also shown in fig- 

ures 4, 5, 6, and 7 . 


Since for a given sideslip angle, the sideslip loading is related 
mainly to the effects of wing deflection and wing angle of attack, both 
of which are linearly related to the airplane normal load factor, (see 
appendix B) the deflection due to sideslip (item ( 7 )) may be expressed 
as 


Z = Z np np 


(5) 


If the effects of aileron, rolling velocity, wing flapping, and 
sideslip are combined, the deflection at any point resulting from the 
unsymmetrical loading only can be expressed by the equation 


Visym = z 5 ail 6 AIL + ^ “Tf + Z np n 3 + V 5 * ^ 


From equation (4) of reference 2, the deflections resulting from 
the symmetrical loadings may similarly be expressed as 


Z sym - Z o + Z n n + Z q3- + Z q^ + Z i 


(7) 


By adding the deflections due to equations ( 6 ) and ( 7 ), the total 
inflight deflections due to both the symmetrical and antisymmetrical 
loadings becomes 


Z = (Z Q + Z ± ) + Z^n + Z^q + Z q q + 
z 6 A il 5ail + Z n Tf n T f + ' L n$ v> $ + 


In the analysis of the roll maneuvers contained herein the unsym- 
metrical part of the deflections was isolated from the total measured 
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deflections by removing the effects of the normal load factor and 
pitching angular acceleration and velocity through the use of all but 
one of the symmetrical deflection coefficients already derived in ref- 
erence 2 . The symmetrical zero-lift and droop deflections Z Q + Z^ 

were not removed from the flight data because the factors which affect 
the zero-lift and droop deflections, such as temperature and possible 
slippages, were suspected to vary from flight to flight.. The deflec- 
tions due to the roll effects plus zero lift and droop were thus 
obtained from the flight data by the equation 

Z'roll = Z - M - Z 44 - Z^q (9) 


During each run the Mach number, dynamic pressure, weight, and 
center-of -gravity position were held effectively constant. Thus, for 
each target in each run the residual deflections, after removing the 
symmetrical portion, may be represented in matrix notation by the 
equation 


{ z 'rollj = ( Z o + Z i) I 1 } + Z 5 a jl { 5 AIl} + 
z nr Ff j^f j + ^ + ^ { p } 

where the columns -j j> on the right-hand side of equation (10) are 

corresponding measured values of &aIl> “Tf* np, and p read from the 

flight records at 0 . 1 -second intervals during the run. The coeffi- 
cients Z Q + Zi, Z 6ail , Z^ , Znp, and Zp for each run were solved 

for by the method of least squares (see appendix A); approximately 
50 data points per run were used. 

The typical time -history plot of the Z' roll deflections given in 

figure 8 shows the 51 points used in the analysis of the roll deflections 
of target number 9 in flight 10 , run 17. Also plotted in figure 8 are 
the 51 data points for the total deflections of target number 9 measured 
during the roll. 

The values of Z Q + Z ± , Zg^, Z^, Z^ , and Zp coefficients 
calculated for each target in each run are presented in tables 3, k, 5, 
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6, and 7 , respectively. Table 8 lists the deflection coefficients for 
each target, the standard error of the coefficients, and the standard 
error of the equation for each target calculated by the method of 
appendix A for one typical run, that is, run 17 of flight 10. In the 
tables no results are shown for optigraph targets numbers 15 and l6 
because the effects of roll on the total deflection of these targets 
were too small to permit a breakdown into the various coefficients. 

Spanwise plots of the wing deflections caused by a unit' aileron 
deflection are shown in figure 9* Similar spanwise plots for the wing 
deflection due to wing flapping, sideslip, and rolling velocity are 
shown in figures 10, 11, and 12, respectively. The variation of these 
wing deflections with dynamic pressure is shown in figure 13 . 

A comparison of the measured wing deflections and those calculated 
from the derived deflection coefficients obtained from the flight data 
is shown in figure 14 for flight 10, run 17 . 


Theoretical Calculations and Comparisons 
With Flight Data 

In the determination of the theoretical curves of wing deflection 
and twist for comparison with experimental results the methods of ref- 
erence 3 were used to calculate the loads acting in each case. The 
lift-curve slopes used in the theoretical calculations were determined 
with data obtained from reference 4. The wing structural stiffness 
distributions were obtained from references 5 and 6. The wing deflec- 
tions resulting from the application of the theoretically calculated 
loads were obtained through the structural influence coefficients of 
reference 1. The theoretical wing twist was calculated by using the 
theoretical structural matrices calculated by the methods of reference 3* 
A detailed outline of the theoretical calculations is given in appendix B. 

Since wing deflections due to unsymmetrical loads are not very large 
as compared to the symmetrical deflections, it was deemed unnecessary to 
obtain experimental and theoretical variations of the various unsymmetri- 
cal deflection coefficients with dynamic pressure and Mach number . As a 
result, the analysis was confined to four roll maneuvers from which only 
one run was selected for comparison with theoretical calculations for the 
same flight condition. 

Comparisons are made in figures 15 to 19 between experimental and 
theoretically calculated deflections and twist due to the aileron loads, 
the wing -flapping inertia loads, the sideslip loads, and the rolling- 
velocity loads on the airplane. The comparisons made pertain to a Mach 
number of O. 7 I, altitude of approximately 35 > 000 feet, and gross weight 
of 122,000 pounds which are the flight conditions of run 17 of flight 10. 
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The comparisons between, calculated and experimental deflections due 
to the aileron axe shown in figure 15- The deflections due to the theo- 
retical aerodynamic load are also shown so as to give an indication of 
the magn itude of the aerodynamic deflections, the inertia deflections, 
and the sum. 

The theoretical and experimental wing deflections due to wing 
flapping are shown in figure l6. The data points shown are the peak 
deflections when the maximum amplitude of the load-factor variation 
caused by wing flapping is unity at the wing -tip accelerometer . . The 
theoretical points, as explained in appendix B, are the double integral 
of the ass ume d accelerations due to wing flapping. Examination of the 
equations of appendix B shows that good agreement between theoretical 
and experimental deflections for those stations inboard of the tip sta- 
tion is a result of the assumption made that the inertia distribution due 
to wing flapping was the same shape as the deflection curve due to wing 
flapping . 

The forces on the airplane resulting from aileron deflection and 
wing flapp ing are generated solely by the wing, but in the case of side- 
slip and rolling velocity there are additional forces on the airplane 
generated by the tail. In the calculation of wing deflections due to 
sideslip and rolling velocity, account was made of the inertia deforma- 
tions caused by the rolling acceleration resulting from the tail loads. 

The experimental and theoretical deflections due to the effects of 
sideslip are given in figure 17 . Shown in figure 17 are the deflections 
caused by the sideslip aerodynamic loads on the wing, the deflections 
caused by the aerodynamic plus inertia loads on the wing, and the total 
wing deflections due to sideslip on the wing and tail. In the theoreti- 
cal calculations the effects of the sideslip loads on the fuselage and 
nacelles were neglected. 

The comparison between the theoretical and experimental wing deflec- 
tions due to rolling velocity is shown in figure 18. In a similar manner 
to that of figure 17 , the theoretical wing deflections due to the rolling- 
velocity wing aerodynamic loads are shown along with the wing deflections 
due to the aerodynamic plus inertia forces in order to give indication as 
to the magnitude of the forces acting. The deflection caused by the 
r oiiing-velocity loads on the tail is also implicitly shown in the figure. 
In the calculation of the loads on the tail (see appendix B), the tail was 
assumed to roll about axis of intersection of the horizontal and vertical 
tail. The calculated effects of the tail loads are considered to be at 
the maximum since no corrections were made for interference effects on the 
tail. The effect of fuselage and nacelle loads on the rolling -velocity 
deflections was also neglected. All of these effects would tend to 
reduce the rolling moment and therefore the wing deflections. 
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The experimental and theoretical wing twist due to the various 
types of wing loadings are shown in figure 19 . The methods by which 
these twists were obtained are outlined in appendix B. 


DISCUSSION 

Deflection Coefficients 


It may be seen from figures 9 to 12 that the deflection coefficients, 
for various wing loads, when plotted against span position for each of the 
spars form well-defined curves. Inasmuch as that the deflection coeffi- 
cients of each target were determined independently of the others, the 
smoothness of the curves indicates that equation (10) which was chosen to 
represent the data probably includes all of the important variables. The 
validity of equation (10) is also substantiated by the agreement shown in 
figure 14 between the measured deflections and the deflections calculated 
using the deflection coefficients and the measured airplane motions. 

The final selection of variables in equation (10) was decided upon 
after first eliminating some of the variables of lesser importance. For 
example, terms expressing deflections due to yawing velocity, sideslip 
effects on the deflected aileron, and the independent effects on the fuse- 
lage and tail were originally included in equation (10). It was found that 
not only was the influence of these components small but also the standard 
errors for these additional coefficients were in some cases as large as 
the coefficients. For the maneuvers considered, the primary variables 
in equation (10) are aileron deflection and wing flapping; sideslip and 
rolling-velocity effects were found to be secondary. 

The Z Q + Z ± coefficients of table 3 obtained from the maneuvers 

analyzed in this report are of the same magnitude as those presented in 
reference 2. The small differences which occur between the two sets of 
data are, as previously mentioned, thought to be attributable to small 
differences in slippages and temperature effects. Theroetical calcula- 
tions, plots of the spar deflections, and variations with dynamic pressure 
are not presented for the symmetrical Z Q + Z^ coefficients because they 

are already covered in reference 2. 

The deflections due to aileron and sideslip (shown in figs. 13(a) 
and 13 (c)) decrease with decreasing dynamic pressure and tend toward a 
value of zero at zero dynamic pressure. This trend is expected since the 
major variable affecting the deflections in each case is the dynamic 
pressure . 


• CONFIDENTIAL 



NACA RM L56C23a 


13 


• • • • • 


• • • • 

• • • 

• • • •• 


• • • 


• • • • • • 

(wwrii®»fAfe»* 




• • 

• • • 


• • 
• • 
• • 


In figure 13(b), the wing -flapping deflections are shown to be inde- 
pendent of dynamic pressure . This fact is true because the inertia forces 
are mainly causing these deflections . 

If the effects of Mach number and flexibility are neglected, the 
wing loads due to unit rolling velocity are theoretically linearly related 
to the square root of dynamic pressure. Thus, the variation with dynamic 
pressure shown in figure 13(d) is somewhat expected. The difference shown 
in the figure between the deflections of a right roll and a left roll is a 
centrifugal effect which is explained by considering that in a left roll 
the left wing deflections due to both the rolling-velocity airloads and 
centrifugal loads are positive and in a right roll the airload deflections 
are negative and the centrifugal deflections are positive. Because of the 
low-slung nacelles and high wing of the B-V7A, calculations indicate that 
for any reasonable vertical location of the axis of roll the overall 
effect of the centrifugal forces is to bend both right and left wing 
upward. 


Theoretical Calculations and Comparisons 

In general, the theoretical methods predicted the experimental 
deflections very well; however, the better agreements between experimen- 
tal data and theory were obtained for the deflections due to the aileron 
deflection and wing flapping which involved only wing forces and which 
were previously mentioned as the most predominate deflections in the roll 
maneuvers. In the case of sideslip and rolling velocity, fuselage and 
tail effects which are difficult to calculate accurately tend to cause 
more disagreement between theoretical calculations and experimental data. 

In the case of the wing deflections due to the aileron, the agree- 
ment between experiment an d theory is good both for vertical deflections 
and wing twist. For wing flapping the vertical deflections are in 
agreement; but the twist is only in fair agreement for the inboard sta- 
tions, and there is a large disagreement at the outboard stations. This 
behavior is believed to be due to the fact that the pitching oscillations 
of the combined outboard nacelle and wing, which are caused by the flapping 
of the wing, were not included in the theoretical twist calculations. 

The deviations of the experimental data from the theoretical sideslip 
deflections, shown in figure 17, are suspected to be mainly a result of 
the faulty assumptions in the theoretical calculations since the experi- 
mental sideslip deflections (table 6) are of good behavior and are consis- 
tently about the same magnitude for all of the runs. As the equations for 
the sideslip effects used (see appendix B) are only approximations to a 
difficult problem, small disagreements can be expected. Also no account 
was taken in the theory for sideslip loads on the fuselage and the nacelles, 
both of which would affect the theoretical answers. In figure 19, the 
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experimental and theoretical twists due to sideslip tend to be in better 
agreement than do the vertical deflections shown in figure 17 . 

The experimental and calculated rolling-velocity deflections shown 
in figure 18 are in fair agreement. The small derivations shown are not 
surprising if the inability to accurately calculate tail loads and the 
neglect of nacelle aerodynamic loads are taken into consideration. The 
scatter in the experimental wing twists due to rolling velocity shown in 
figure 19 (c) axe probably a result of the small rolling deflections 
obtained in the flight maneuvers (that is, on the order of 1 in. at the 
wing tip ) . 


CONCLUDING REMARKS 


The inflight wing deformations on the B-kjA airplane have been ana- 
lyzed for the effects of the various types of wing loadings and the 
results are presented as deflection coefficients. Theoretical results 
combining existing methods of aerodynamics and stress analysis are also 
presented and show good agreement with the experimental data. Because 
of the particular aerodynamic and inertia characteristics of the Boeing 
B_k7A wing, the wing deflections due to roll effects were relatively 
smaller than the symmetrical deflections previously reported in NACA 
RM L5 1 +K24a . Even though these deflections are small, a detailed analysis 
was possible due to the high measuring accuracy of the optigraph system. 


Langley Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 
Langley Field, Va., March 13, 1956. 
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APPLICATION OF LEAST -SQUARES METHOD TO EVALUATING 
DEFLECTION COEFFICIENTS 


This appendix outlines the method used to reduce the wing deflections 
of this report to coefficient form. Since the method has a wider applica- 
tion, it is also illustrated in detail by obtaining numerical results for 
one of the time histories given in the body of the report. 


The equation to which the least-squares methods are applied is equa- 
tion (10) of the text, namely, 

= (Z 0 + z i) {l} + | 5 AIl|’ + 

Zn T f {^f } +Z np{ np } + ^j 1 } 

are corresponding values of 2. ! ro n> b ATT > 

nrp_^ np, and p obtained from the flight records at various times 

during the run. Example values obtained from figure 5 and equation (9) 
for these columns are listed in the following table : 


{z'roll} 


where the columns 


•{ 


Time, 

sec 

5 ail> deg 

“Tf 

np, deg 

„ radian 
P> sec 

Z 'roll ( f: pom 

eq. (9)) 

0 

-0.082 

0.005 

0.57^ 

-0.004 

23.144 

1.0 

-9.316 

-.254 

.658 

-.146 

19.719 

2.0 

-9.481 

-.462 

-.624 

-.297 

19.910 

3.0 

-6.986 

-.016 

-2.302 

-.273 

19.604 

4.0 

0 

.517 

-3.136 

.038 

21.274 

5-0 

-.415 

.308 

-.382 

.082 

21.818 


When there are more equations than unknowns the equations must first 
be normalized in accordance with least-squares procedures given in refer- 
ence 7 . ' The normal equations from which a solution for the unknown 
deflection coefficients were obtained aa?e as follows : 
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< z ° + Zi)H + z 6ail £ »ail + % f £ “Tf + 


Znp ^ n^ + Zp ^ P = 

y Z'roll 


(Al) 

( z o + z i) ^ 5 AIL + Z 6 ML £ 

Sail 2 + S f X 

“Tf 6 AIL + 


^3 ^ n P 5 AIL + ^ ^ p5 AIL 

= y Z' r oll 6 AIL 


(A2) 

(Z 0 + Zi) Y °T f + Z 6 ajl ^ 

& AIL n T f + ^ 

Z “Tf 2 + 


Znp . ^ ^"Tf + h ^ pn T f = 

; y Z ' r oll n T f . 


(A3) 

(Zo + z i) y np + Zq^ y 

8 AXL n P + Zruj,^ y 

' nr f np + 
1 


Znp ^(np) 2 + pnp = ^ 

) z ' roll np 


(Ak) 

(z D + ^ y p + z 6aiJj y 

j 6 ailp + z * If y 

1 

m-jP + 


Znp y npp + Zp y P 2 = 

y Z'rollP 


(A5) 

The value of N in the preceding equations is the number of 

sets 


of data available; in this example, six sets of data are available. In 


these equations the quantity 


^ 5 AIL is simply tlie summation of the 


5 att column in the preceding table. The quantity 


y 8 ail 2 is sim P 1 y 
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the summation of the square of each number in the 5 ATT column. The 
quantity ^ 5 A - rTi n Tf is obtained by multiplying, at each listed time, 

the value on the bpjj, column by the corresponding value in the 
nrji^ column and summing all of the obtained products. A similar proce- 
dure is followed in determining all of the other summations in the pre- 
ceding normal equations. 

After the normal equations have been determined, the solution for 
the unknown coefficients Z Q + Z±, ^ ah/ ^T ’ ^P ma ^ r 

obtained. The normal equations for the example are 

6(z 0 + Zi) - 26.280 Z 6ail + 0.098z rlTf - 

5.212Z n p - 0.592Zp = 125.469 (a6) 


-26.280(z o + Zj.) + 255-660362 Z 5ail + 6.730032Z nf ^ + 

15 • 979450Z n p + 6 . 0488i3Zp = -520.374736 (A7) 

0.098(z o + Zi) + 6.730032 Z 6ail + 0.640394z nT ^ - 

1.578llOZ np + 0.223588Zp = 3.313612 (a8) 

-5.212(Z 0 + Zi) + 15.979450 Z &ail - i^noz^ + 

16.431440ZJJP + 0.569510^ = -106.342230 (A9) 

-0.592(Z o + Z ± ) + 6.048813Z 5ail + 0.223588z n ^ + 

0.5695l0Z np + 0.l92238Zp = -11.454072 (AlO) 
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The preceding normalized equations were obtained from only six sets 
of data and are presented here only to demonstrate the method; actually 
the coefficients were determined from 51 sets of data in order to obtain 
more accurate values of the deflection coefficients. Further discussion 
is based on the normal equations obtained from 51 sets of data points for 
flight 10, run 17 of figure 5 as follows: The normal equations in this 

case were 


5l(z 0 + Zi) - 233.831 Z 5ail + 3.066z_ - 

55.327^ - 5.803Zp = 1057.730 (All) 

-233. 831 (Z Q + Zi) + 2092. 57422 3 Z Qail - 11.527022^ + 
ll6.80782iz np + 54.0l7682Zp = -4437.738398 (A12) 

3.066(z o + z ± ) - 11.527022Z &AIL + 7.743859Z nTf - 
8.i39764z np - 0.447479Zp = 47.507306 (A13) 

-55-327(Z 0 + z i) + h6.807821Z 5ail - 8.l39764z rlTf + 
l6l.25l305Znp + 5.324052Zp = -1185. 655375 (Al4) 

-5.803(z o + Zi) + 54.017682 Z &ail - 0.447479^ + 
5.324052Z np + 1.749 17 3Zp = -109.307618 (A15) 


which when solved simultaneously give 


(z Q + Zi) = 23.047 z 6ail = 0.494 Z- = -2.218 
z np = 0.163 Zp = -2.346 
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Substitution of the coefficients in equation (10) gives: 

Z' ,, = 23.047 + 0.494Z r - 2.2l8n rp + 0.l63np - 2.346p (Al6) 

roil ' Sail 


The solution to the normal equations presented in this paper were 
obtained through the use of Crout's method (ref. 8), because it was 
easily applicable to computing -machine use. 


An estimate of the ability of equation (Al6) to fit the 
deflections is given by 


Standard error of fit = 



measured 


(A17) 


where N equals the number of equations (51 in the case cited) and U 
eq ual s the number of unknowns (5). 


The summation of deviations squared term, 

from : 


^ (dev^), was obtained 


dev 2 = ^ ( z 'roll ) 2 " ( Z o + Y z 'roll “ 

z &AIl Yj Z ' roll5AIL ” Zn T f Y Z ' rolinT f ' 

Znp Y Z ’roll n P “ h Y Z ' rollP (Al8) 


These s umraations are given for N = 31 in equation (All) to ( A15 ) , 
with the exception of z ' roll 2 which were obtained by summing the squares 

.of the 51 values of Z'roll* The standard - error of fit for target 9 was 

+O.I76 inch, which indicates that the average deviation of the measured 
deflection data points from the determined equation is approximately 
+O.I76 inch. 

The standard error of each of the coefficients was obtained by the 
equation 

jstandard error of coefficients | = (Standard error of equation) 

(A19) 
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The values of the /T coefficients were determined from the numbers 
in the matrix of the normal equation (All) to (A15), as follows. 

The normal equations (in matrix form) for the 51-P°i n t example pre- 
sented previously are: 


51 

-233.831 

3.066 

-55.327 

-5.803 


(Zo + Zi? 


1057.730 

-235.831 

2092.574223 

-11.527022 

116.807821 

54.017682 


Z6 ail 


- 4437.738398 

3.066 

-11.527022 

7.743859 

-8.139764 

^0.447479 

< 

Vrf 

> = < 

47.507306 L 

-55.327 

116.807821 

-8.139764 

161.251305 

5 . 324052 




-1185.655375 

-5.803 

54.017682 

-0.447479 

5,324052 

1.749173 


. ^ > 


-109.307618 


(A20) 


or 


M 


z Q + z± 


Y Z ’roll 

Z &AIL 


^ z 'rOll 5 AlL 

< Zn Tf 

► = < 

^ z 'roll n T f > 

Z nf3 


^ z ' roll n P 

v „ 


Z 'rollP 

V J 


(A21) 


The inverse of matrix M is (giving only the diagonal elements) 
lo. 0910146 — — 


[A ]" 1 - 


0.0048902 


O.I 36797 O 


0.0164597 


3.7846045 


(A22) 
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or 



C 22 



— — 


(A23) 


c 


55 


Thus the standard error S of the coefficients denoted by the subscripts 
for this example were 


f y 

s (z 0 + z i) 


'fiTi 


±0.053 

Zb ail 


J C 22 


±0.011 

4 S Zn Tf 

>■ ~ (^equation ) A 

/ c 55 

► = < 

±0.065 > 

Sz n P 


\fikk 


±0.022 

% 

'k J 


|/^55 

V. J 


±0.3^2 

t 4 


(A2M 
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In the theoretical calculations presented in the body of the report 
the effects of the wing airload, the effects of the rolling acceleration 
due to the wing airloads, the fuselage airloads, and the airplane tail 
loads have to be considered. The general procedure was to calculate the 
wing airloads and the rolling acceleration resulting from these loads by 
the methods of reference 5 and, when necessary, to approximate the rolling 
acceleration due to the fuselage and tail loads by the methods included 
in this appendix. 


S1MB0IS PERTINENT TO APPENDIX B 


A' ' amplitude of oscillation, in. 

0 m aileron-section pitching moment due to unit aileron deflection 

& AIL 

g acceleration due to gravity, 386-4 in. /sec 2 

I airplane rolling moment of inertia, 490 X 10^ lb-in. ^ (obtained 

from ref. 9 ) 

K c aileron-section effectiveness factor 

°AIL 

m 0 wing-section lift-curve slope 

r spanwise distance to vortex center line, in. 

Q 

S total wing area, in. 

V true free-stream velocity, in. /sec 

a wing-section streamwise angle of attack, positive leading 

edge up 

t time, sec 


dn 


wing-section total angle of attack per unit load factor 
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d0 Y 

— — slope of elastic axis deflection curve per unit load factor 

dn 

A 


0 ) 


wing sweep angle , positive for sweepback, deg 
circular frequency of the wing flapping 


Matrix notation: 


M 

LKgJ 

{'} 

N 


structural matrix expressing the wing angle of twist caused 
by the nacelle pitching moment (ref. 3) 

geometry matrix to determine rolling moment associated with 
known wing loads 

wing-section lift, lb/in. 


angle of attack induced on the wing due to the fuselage 
overvelocities caused by sideslip 


[Si + S ± ] aerodynamic downwash matrix calculated from the geometry of 
the wing and fuselage vortices 


M 

ft 


structural matrix relating wing twist to the wing loads 


diagonal matrix in which diagonal elements are equal to unity 


/ 




wing-section streamwise angle of attack in the flight condi- 
tion due to all effects other than the angle of twist due 
to 0 . 25 c airloads 


•ta,. 1 wing-section angle of attack caused by wing twist due to aileron 

y ‘"Q flTT. j pitching moment caused by unit aileron deflection 


H 


{-} 


wing angle of attack caused by the rolling acceleration inertia 
deformations 

structural matrix expressing the wing deflections due to the 
wing -quarter -chord loads 

column matrix expressing the centrifugal deflection of the 
wing caused by a unit rolling velocity at n = 1 
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structural matrix expressing target deflections per unit 
rolling acceleration 

matrix of wing target deflections 


matrix of wing -flapping normal load factor at 
targets, when = 1 


the op ti graph 


q* dynamic pressure 

The theoretical unsymmetrical aerodynamic loads are given by equa- 
tion (24) of reference 3 as 


[l]- q" S W][^][ Sl + S i] ' [M {*} ' {°^} (B1) 


In order to account for rolling acceleration flexibility effects, 
where the rolling acceleration is a function of the airloads on the 
wing, the equation 


|_ k gJ {* } = x xP 


(B2) 


was solved simultaneously with equation (Bl) so as to give the wing 
loads and the rolling acceleration associated with each flight condition. 


The effective wing angles of attack 


w 


due to aileron, side- 


slip, and rolling velocity were determined separately and inserted in 
equation (Bl) so that each solution gave the wing airload and the rolling 
acceleration resulting from that angle -of -attack distribution. 


The wing bending and torsional stiffness distribution of references 5 
and 6 were used in the determination of the structural matrices in equa- 
tion (Bl). The wing was treated as a cantilever beam fixed at the inter- 
section of the wing 38 -percent-chord line and the fuselage wall. The 
geometry and downwash matrices of the above equation were calculated for a 
vortex system of 10 vortices per semispan. The spanwise location of these 
vortices may be deduced from the twist comparison figure (fig. 19) where 
the theoretical twist data points are shown. The values of the section- 
lift-curve slopes mo were determined through the use of the above equa- 
tion as applied to the wind-tunnel symmetrical data of reference 4 and 
includes correction to the Mach number by use of the Prandtl-Glauert 
equation. 
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In the determination of the aileron airload on the flexible wing 


j a gj = | K 5 ai L | 6 AIL + p AIL + 


:m& AIL 


3 AIL 


The aileron-section effectiveness factors. 


i 


| k& ailj 


(B5) 


were obtained 


was determined by using the 


/ OjO 

obtained from reference 6 in conjunction with the 


from reference 10. The matrix 
values, of 

proper structural matrix determined by the method previously mentioned. 
The ct£ column matrix was calculated by applying the rolling accelera- 
tion inertia distribution to the wing through theoretical structural 
matrices . 

In the case of rolling velocity p, the equation 


{“g| * {v} p+ 




again includes the effects of inertia on the flexible wing. 
The sideslip a CT ^ distribution was calculated as 


M ■ {& 


■M 


.p tan A + 


{ S °h} p + {^} pcosA+ 


(B5) 


which is the result of the flexible swept wing and fuselage sideslipping 
at n = 1 and which includes the effects of wing dihedral and fuselage 
overvelocities calculated by the methods of references 11, 12, and 13 . 

In the case of sideslip and rolling velocity, where the associated 
rolling acceleration of the airplane is affected by the tail loads, the 
wing deflection caused by the tail loads was handled separately and 
superimposed on wing effects to obtain the wing deflection due to the 
combined effects of wing and tail. The vertical -tail-load effect in 
sideslip was calculated by using a C L(j of 0.045 (estimated from pre- 
liminary strain-gage data obtained in the Flight Research Division of 
the Langley Aeronautical Laboratory) and by assuming the rolling axis 
to be through the airplane center of gravity and parallel to the fuse- 
lage axis. The horizontal tail was assumed to act in the same manner 
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as the wing in sideslip. The angle of attack of the horizontal tail was 
calculated from the airplane tail-off pitching moment at n = 1. 

In calculating the effect of tail load on rolling velocity deflec- 
tions the tail was assumed to be rolling about an axis passing through 
the intersection of the vertical and horizontal tail. The wing c^ dis- 
tribution due to rolling velocity was assumed to act on the tail. 

Once the total airloads and rolling accelerations were determined 
for the flight condition, the theoretical-angle -of -attack change due to 
wing twist was calculated from the equation 

<JaJ^ = [s 2 ]|z| + |p (B6) 

In the case of wing twist due to aileron, an additional term 

Ja P 6 att w as included in equation (b 6) to account for wing twist 

| m SAlLj ■ 

due to aileron pitching moment. 

In the calculation of the theoretical angle of twist due to wing 
flapping, an inertia distribution which had a value of unit load factor 
at the wing tip and was the shape of the wing -flapping deflection curve 
was applied to the wing through the use of a proper structural matrix 
calculated by the method previously mentioned. No account was taken of 
pitching oscillations of the wing or nacelle caused by wing flapping. 

The deflections were determined from the equation 


{4 - ini {»}■ 

where the matrices ||^z|| 
ence coefficients of reference 1. 




(B7) 


were determined from the influ- 


The wing deflection due to aileron was calculated by using equa- 
tion (B7) and altering the || 5 z|| niatrix to account for the shift in 

position of the lifting line from the quarter-chord line at the aileron 
in accordance with reference Ik . A term w p2 was added to equa- 
tion (b 4) in the calculation of rolling -velocity wing deflections, which 
accounted for the centrifugal effects on the deflected wing. In the 
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matrix, the vertical location of the air- 


plane center of gravity was assumed to he at Z M = 56 inches, obtained 
from reference 9* 


The deflections associated with wing flapping were theoretically 
checked by assuming sinusoidal oscillations and calculating the ampli- 
tude of oscillation from the equation 


— z = A 1 ^ sin o)t 
dt 2 


nfg 


o 

A'orsin cut 


where o> was measured from the time-history plots of n T^* ^ ie 

wise shape of the wing acceleration distribution due to wing flapping n^ 

was assumed to be of the same shape as the wing deflection curve associ- 
ated with rif • 

Thus, at the maximum amplitude where sin cot = 1 

{4 
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TABLE 1.- FLIGHT CONDITIONS 


Flight 

Run 

Maneuver 

Mach 

number 

Dynamic 
pressure, 
lb/sq ft 

Altitude, 

ft 

Aircraft 

weight, 

lb 

Position 
of c.g., 
percent 
M.A.C . 

8 

15 

Left roll 

0.60 

125 

34,500 

116,000 

16. 04 

10 

17 

Left roll 

• 71 

177 

34,900 

122,000 

22.48 

10 

18 

Right roll 

• 71 

171 

35,300 

121,000 

25.79 

12 

14 

Right roll 

.6o 

130 

34,800 

117,000 

13.50 
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a Read manually as slope of velocity trace. 
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TABLE 3-- ZERO-LIFT PLUS DROOP WING DEFLECTION COEFFICIENT Z Q + Z ± 


Flight 

Run 

Zero-lift plus droop wing deflection coefficient, in., at target - 

9 

10 

11 

12 

15 

14 



■8 

15 

21.897 

22.727 

13.100 

13.009 

7.5OO 

8.020 

10 

17 

23.047 

24.464 

13-516 

14.388 

7.387 

8.560 

10 

18 

23-559 

24.695 

13.464 

n.517 

7.501 

8.623 

12 

14 

19.205 

20.544 

11.986 

12.283 

7.130 

7.551 


TABLE 4.- AILERON WING DEFLECTION COEFFICIENT Z &AJL 


Flight 

Run 

Aileron wing deflection coefficient, in. /deg, at target - 

9 

10 

11 

12 

13 

14 

8 

15 

O.368 

0.492 

0.179 

o'. 209 

0.054 . 

0.093 

10 

17 

.494 

.578 

.237 

.296 

.096 

.143 

10 

18 

.424 

.528 

.206 

.255 

.076 

.121 

12 

14 

.308 

• 559 

.145 

.181 

.059 

.087 


TABLE 5 . - WING-FLAPPING WING DEFLECTION COEFFICIENT 


Flight 

Run 

Wing -flapping wing deflection coefficient, in./nTf, at target - 

9 

10 

11 

12 

• 13 

14 

8 

15 

-2.219 

-2.136 

-1.176 

-1.241 • 

-0.493 

-0.6l4 

10 

17 

-2.218 

-2.247 

-1.190 

-1.3<*9 

-.510 

-.695 

10 

18 

-2.131 

-2.225 

-1.113 

-1.196 

-.461 

-.629 

12 

14 

-2.310 | 

-2.194 

-1.269 

-1.342 

-.560 

-.667 


TABLE 6.- SIDESLIP WING DEFLECTION COEFFICIENT 


Flight 

Run 

Sideslip wing deflection coefficient, in. /deg, at target - 

9 

10 

11 

12 

13 

14 



8 

15 

0.111 

0.128 

0.074 

0.074 

0.035 

0.049 

10 

17 

.163 

.192 

.111 

.122 

.048 

.029 

10 

18 

.152 

*157 

.094 

.094 

.037 

.042 

12 

14 

.126 

• 073 

.065 

.076 

.009 

.031 


TABLE 7. - ROLLING -VELOCITY WING DEFLECTION COEFFICIENT 7 ^ 


Flight 

Run 

Rolling -velocity wing deflection coefficient, in. /radian/ sec, 

at target - 

9 

10 

11 

12 

13 

14 

8 

15 

-2.540 

-5.252 

. -1.278 

-1.482 

-0.387 

-0.999 

10 

17 1 

-2.346 

-2.674 

-1.255 

-1.505 

-.788 

-.753 

10 

18 

-1.411 

-2.708 

-.872 

-1.259 

-.317 

-.633 

12 

14 

-1.283 

-1.777 

-.797 

-1.102 

-.524 

-.681 
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Figure 1. - Photograph of test airplane 
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Figure 5 .- Time-history data for flight 10, run IT, left roll, Mach num- 
ber 0 . 71 , altitude 3^,900 feet. 
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Figure 6.- Time-history data for flight 10, run 18, right roll, Mach num- 
ber 0.71, altitude 35 >300 feet. 
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Airplane normal load factor, n p, radians per second _ 

Wing-tip normal load factor, nj P> radians per second per second ^9 * 



Mach num' 


Wing-flapping load factor, n T q, radians per second. P* degrees; degrees 

f radians per second per second 
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o Data point 


0 .4 .8 1.2 1.6 2.0 2.4 2.8 3.2 3.6 4.0 4.4 4.8 5.2 

Time, seconds 

Figure 8. - Typical time history of target 9 deflection data, before and 
after symmetrical deflections were removed, showing 51 data points 
used in analysis. 
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Figure 9* - Flight-test front-spar deflections due to aileron deflection 
of 1°, altitude approximately 35>000 feet. (See table 1. ) 
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Figure 10.- Flight-test front-spar deflections due to wing-flapping load 
factor, nT f = 1, altitude approximately 35>000 feet. (See table 1. ) 
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Figure 11.- Flight-test front-spar sideslip deflection per np, altitude 
approximately 35*000 feet. (See table 1. ) 
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Figure 12.- Flight-test front-spar deflection due to rolling velocity of 
one radian per second, altitude approximately 35*000 feet. (See 
table 1. ) 


CONFIDENTIAL 














k6 


CONFIDENTIAL* ** 


NACA RM L56C23a 


Q) 

a> 

5m 


SO 


CD 

r r3 

5m 

CD. 

P- 


CQ 

CD 

X 

O 

C 


I— I 


CO 

CS3 


2.8 
2.4 
2.0 
1.6 
1.2 
.8 
.4 

° 0 .1 -2 .3 .4 .5 .6 .7 .8 .9 1.0 

Y’/b'/2 








— n 

Theory, 



— Fro 
ReJ 

nt s.] 
a.r sp 

Ex] 

Dar 

ar 

Derim 

c 

u 

ental 

ai 

rload + ii 

V 

V 

1 

lerti 

A 

V 


t 

Theor 

•y,ai: 

rload 

onlj 


/ 

/ 


r \ 

>" 













% 


. - 


^ 


f 


§ 




Figure 15. - Comparison of experimental and theoretical aileron wing 

deflection coefficients. 
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Figure l6.- Comparison of experimental and theoretical wing-flapping 

deflection coefficients. 
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Figure 17.- Comparison of experimental and theoretical sideslip wing 

deflection coefficients. 
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Figure 18. - Comparison of experimental and theoretical rolling-velocity 

wing deflection coefficients. 
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Figure 19*- Comparison of experimental and theoretical streamwise wing 

twist. 
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